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The constraints on the asymptotic behavior of the conformal factor and conformal extrinsic curva- 
ture imposed by the initial value equations of general relativity on constant mean extrinsic curvature 
(CMC) hypersurfaces are analyzed in detail. We derive explicit formulas for the Bondi-Sachs energy 
and momentum in terms of coefficients of asymptotic expansions on CMC hypersurfaces near future 
null infinity. Precise numerical results for the Bondi-Sachs energy, momentum, and angular momen- 
tum are used to interpret physically Bowen-York initial data on conformally flat CMC hypersurfaces 
(N- similar to that calculated earlier by Buchman et al. 
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I. INTRODUCTION 



\^ • Recently, there has been increased interest in formulations for numerical relativity based on conformal compact- 
ification [2| in which the calculational grid extends to future null infinity {J^'^), where the gravitational radiation 
amplitude can be read off unambiguously, with at most numerical errors, and where no dynamical boundary condi- 
^ ■ tions are necessary. In principle this can be done with Cauchy characteristic matching methods, but these have not 
been implemented successfully with non-trivial dynamics. Cauchy characteristic extraction methods can extrapolate 
^ • from the outer boundary of a conventional Cauchy code to determine waveforms at , but do not eliminate errors 
{^J^ in the Cauchy development deriving from inexact boundary conditions at a finite radius. The characteristic methods 
I— are reviewed by J. Winicour Q. Our focus here is the initial value problem on hyperboloidal spacelike hypersurfaces, 
^ [ and specifically the case of constant mean extrinsic curvature (CMC) hypersurfaces H-Q- The vanishing of the con- 
J> ■ formal factor Q at ,y~^ accounts for the singular behavior of the physical spacetime metric at J^~^ in compactificd 
' coordinates. The conformal geometry, in which is an ingoing null hypcrsurfacc, is regular in a neighborhood of 
^\ , , and a CMC hypersurface of the physical spacetime is spacelike in the conformal spacetime out to and including 

lO ' its intersection with a 2-surface with spherical topology which we denote by Regularity conditions relating 
. the 2D extrinsic curvature of ^+ as embedded in the CMC hypersurface to the conformal extrinsic curvature need to 
be satisfied at but when imposed in the initial data are automatically preserved by the evolution equations d, Q. 
The constraint equations, given suitable gauge conditions, determine the leading behavior of the conformal factor and 
the conformal extrinsic curvature of the hypersurface in the neighborhood of .y'^ in terms of asymptotic gravitational 
wave amplitudes @. If the physical mean extrinsic curvature is not too large, the dynamics of the sources (black holes, 
neutron stars) takes place where the CMC hypersurfaces are not that different from the hypersurfaces of conventional 
3-1-1 methods. 

Conformally flat initial data on CMC hypersurfaces were considered in Ref. [l| . The conformal momentum constraint 
cd [ equations have the same form as they do on maximal hypersurfaces, and with conformal flatness admit a class of 
analytic solutions which are slight generalizations of the well-known Bowen-York solutions [l^l often used for single or 
multiple black hole initial data on maximal hypersurfaces. The Hamiltonian constraint provides an elliptic equation 
for the conformal factor which is degenerate at J^~^. Despite this degeneracy, Ref. [l| obtained numerical solutions 
without much difficulty using the standard spectral elliptic solver of the Caltech-Cornell-CITA SpEC code [ll|, [l^ . 
The degeneracy constrains the leading terms in the expansion of fl about 

An essential part of the physical interpretation of these solutions is knowing precisely the total energy, linear 
momentum, and angular momentum of the system as coded in the asymptotic behavior of the spacetime metric at 
The standard Arnowitt-Descr-Misner (ADM) formulas [l^ for these quantities only apply on asymptotically 
fiat slices at spatial infinity. The Bondi-Sachs energy and momentum at J^~^ [13, [3 relevant quantities for 

CMC hypersurfaces. Ref. [l| did not fully address this issue, making only some rather crude estimates of the total 
energy, momentum, and angular momentum with limited validity and for the most part with uncertain errors. 

While there is an extensive literature dealing with the problem of extracting these global physical quantities near 



or at at future null infinity (see, e.g., the review [16[), part of which deals specifically with CMC hypersurfaces [17|, 
we see practical difficulties in implementing many of these procedures. Some of them require the choice of a reference 
spatial metric and extrinsic curvature. Determining the reference quantities with the precision necessary to obtain 
unambiguous results requires considerable effort, in general. Formulas based on the asymptotic behavior of the 
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Weyl tensor (e.g., can be useful computationally, but an analysis based on geometrically defined spatial 

coordinates gives more insight into the asymptotic geometry of CMC hypersurfaces and its relationship to the Bondi 
expansion on null hypersurfaces. 

We first consider the problem of calculating the Bondi energy and momentum in general asymptotically flat space- 
times foliated by CMC hypersurfaces. Following Appendix A of Ref. Q, we analyze the asymptotic geometry using 
Gaussian normal spatial coordinates tied to the J'^ 2-surface. This gauge choice allows a relatively simple charac- 
terization of the asymptotic behavior of the spacetime metric and extrinsic curvature on a CMC hypersurface based 
on asymptotic solutions of the constraint equations. Constructing the asymptotic coordinate transformation from the 
CMC-based coordinates to Bondi-Sachs null coordinates gives us the Bondi-Sachs mass aspect in terms of geometric 
objects defined on a single CMC hypersurface. The monopole and dipole moments of the Bondi-Sachs mass aspect 
are the Bondi-Sachs energy and linear momentum, respectively. In the simple case of conformally flat initial data, 
the mass aspect is just a sum of contributions from the asymptotic expansions of the conformal factor and conformal 
extrinsic curvature. 

The Bowen-York solutions for the conformal extrinsic curvature contain parameters such as boost and spin vectors, 
which on maximal hypersurfaces are simply related to the ADM momentum and angular momentum. On CMC 
hypersurfaces, the relationship of the Bowen-York parameters to physical momenta and energies is more complicated, 
but we are able to obtain analytic expressions for the direct contributions of the conformal extrinsic curvature to the 
Bondi-Sachs energy and momentum in terms of the Bowen-York parameters. Unlike on maximal hypersurfaces, on 
conformally flat CMC hypersurfaces the coordinate displacement of the black hole from the center of the coordinate 
sphere representing enters in a non-trivial way. The conformal factor contributions to the Bondi-Sachs energy 
and momentum are extracted from the numerical solution of the Hamiltonian constraint equation. We also derive an 
analytic expression for the total angular momentum of the system which only depends on the Bowen-York parameters, 
without any contribution from the conformal factor. We discuss the physical interpretation of some representative 
examples of Bowen-York initial data on CMC hypersurfaces similar to those of Ref. [l| , and show how to construct 
initial data approximating a circular-orbit binary black hole system. 

The Wald sign convention for the extrinsic curvature is adopted, so that the mean extrinsic curvature K is positive 
for hypersurfaces with diverging future-directed normals, the hypersurfaces extending to future null infinity rather 
than to spatial infinity (K = 0) or to past null infinity [K < 0). Our notation is an amalgam of that of Refs. P, 
and very similar to that in Appendix A of Ref. @, but be alert to occasional deviations. We assume vacuum in a 
neighborhood of future null infinity. The discussion of the asymptotic geometry on general CMC hypersurfaces is 
in Sec. |IT1 In Sec. Illli we use the lapse and shift which preserve the CMC hypersurface condition and the Gaussian 
normal spatial coordinate condition to facilitate the asymptotic transformation to Bondi-Sachs coordinates. The final 
expression for the mass aspect depends only on the geometric properties of a single CMC hypersurface. Numerical 
results for the generalized Bowen-York initial data on conformally flat CMC hypersurfaces are presented in Sec. IIVI 
The implications of our results for astrophysically interesting initial data are summarized in Sec. El In Appendix [XI 
we derive the analytic expressions for the conformal extrinsic curvature contribution to Bondi-Sachs four-momentum 
and in Appendix [B] the analytic expression for the total angular momentum. 

II. ASYMPTOTIC BEHAVIOR ON CMC HYPERSURFACES 

The fact that the Hamiltonian constraint equation for the conformal factor and the momentum constraint equation 
for the extrinsic curvature are degenerate at J"^ means that the first few terms in a power series expansion of these 
quantities away from ^+ are uniquely determined from the expansion of the conformal spatial metric, in spite of the 
elliptic character of the equations and therefore the generally nonlocal character of the solutions. In our derivation 
of formulas for the Bondi-Sachs energy-momentum 4- vector, we will not assume the conformal spatial metric is flat, 
though all of our numerical examples are for the conformally fiat case. The expansion parameter defined in Ref. [l[ is 
only meaningful for a fiat conformal spatial metric, so the analysis here follows the derivation of asymptotic behavior 
presented in Appendix A of Ref. Q , which in turn rather closely follows that of Ref. @ . It is based on a particular 
choice of spatial coordinates, Gaussian normal coordinates, in which the "radial" coordinate z is the conformal 
proper distance inward from ^+ along normal spatial geodesies in the CMC hypersurface considered as a slice of the 
conformal spacetime. As in Ref. Q, we assume that the intrinsic geometry of ^+ is that of a 2-sphcrc of fixed area in 
the conformal geometry, which can be enforced in any asymptotically flat spacetime by suitable boundary conditions 
on the initial conformal gauge and its subsequent evolution, at least if the gauge evolution variables are determined 
by elliptic equations, as in Refs. Angular coordinates are propagated along the normal geodesies starting 

from standard polar angles on the two-sphere. While the actual computational coordinates will be different, in 
general, the coefficients in the expansions are expressed in terms of quantities defined covariantly and therefore can 
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be computed in any coordinate system. An equality which holds only at is indicated by a dot above the equal 
sign. 

We just give the key definitions and results here, and refer the reader to Rcf. Q for details of the derivations. The 
conformal spacetime metric g^j/ is related to the physical spacetime metric g^^ by 

and a tilde will be used generally to distinguish quantities associated with the conformal, rather than the physical, 
geometry. The conformal factor is required to vanish on where it cancels the infinities in the physical metric 
due to the compactified spatial coordinates and the CMC hypersurfaces becoming asymptotically null, and is strictly 
positive in the interior. The general form of the conformal spatial metric in Gaussian normal coordinates is 

cjijdx^dx^ = dz^ + JiABdx^dx^ , (1) 

with z = at J^+. A metric function ^ is defined such that det Iiab = S,~'^ sin^ 9. At , Kab = S.o'^hAB, with Hab 
the standard metric of the unit two-sphere in polar coordinates. We require that be independent of angle. Then 
the 2D scalar curvature of is 2Cq. The only way to preserve this property of ^+ is to demand that 9(^0 = 
during the evolution as a condition on the evolution of the conformal gauge. We can keep the CMC time coordinate 
t equal to a Minkowski retarded time coordinate u at J^"*" by imposing the boundary condition a == K/(3^o) = ao on 
the elliptic equation for the conformal lapse which follows from dtK = 0. 

A general form for the expansion of Hab in powers of z away from is then 



hAB = £, 



hAB - 2xABZ + (x^'^XcohAB - ^ab) + 0{z'^) . (2) 



Here xab and ipAB are traceless, symmetric tensors on the unit sphere with indices lowered and raised by Hab and 
its inverse. Note the identity that Xa Xcb ^ ^x'^^XcohAB- As argued in Ref. we believe it is consistent and 
physically appropriate in the context of gravitational radiation generated by the internal dynamics of isolated systems 
not to allow any polyhomogeneous termSj^ that is, terms containing powers of logz as well as powers of z, to this 
order in the expansion. The expansion of Hab is closely related to the expansion of the angular part of the metric in 
Bondi-Sachs coordinates. Ref. [20| showed that the coefficient of the leading polyhomogeneous term of order logr 
or r^^ log r in the Bondi expansion is a constant of the motion, so if these terms are absent initially they never appear. 
If these polyhomogeous terms are not present in the Bondi expansion, the corresponding terms are not present in 
Eq. [2] The inverse of the angular part of the metric is 



^AB ^ ^2 ^^AB ^ 2^AB^ ^ [x^^'xCDh^^ + j Z' + 0(Z^)J . (3) 

The angular metric functions can be related to the extrinsic curvature of the constant-z two-surfaces as embedded 
in the 3D conformal geometry of the CMC hypersurface. The trace of the 2D extrinsic curvature is k = 2dz (log C) , 
and the traceless part is 

Ki^xi + ^4z + oiz^). 

The expansion of as a power series in z can be expressed in terms of the expansion of k, 

K — ko + kiz + K2Z'^ + O(z^), 



with the result 



logC ^ log Co + ^koz + \kiz^ + ^^22^ + 0{z''). (4) 



Satisfying our condition on the intrinsic geometry of J"^ requires that be a constant, independent of both time and 
angle, but all of the coefficients in the expansion of k will in general depend on the time t and the angular coordinates 
x"^ . The 2D extrinsic curvature is a tensor which can be calculated in any coordinate system by solving the geodesic 
equations in the conformal spatial geometry. 

A straightforward calculation gives the conformal 3D Ricci tensor components. The result is Eq. (A6) of Ref. [1], 

R\ = dzk - k-ABk^^ - , 

R^A 2^\A - fi'^ A\c, (5) 
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with 
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Up, 



so 



Kq 



X iCD 



(6) 



The I symbol denotes a covariant derivative with respect to the metric Hab of the constant-z two-surface, while the ' 
symbol denotes a covariant derivative on the unit sphere. The curvature scalar evaluated at is 
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R = 2^0^ + 2^1 - 



X XAB- 



(7) 



We define the 3D conformal extrinsic curvature Kij to have the same relationship to the time derivative of the 
conformal spatial metric, the conformal lapse, and the shift as the physical extrinsic curvature Kij has to the time 
derivative of the physical spatial metric, the physical lapse, and the shift. This corresponds to the convention in 
Rcf. and the usual convention in the literature, but differs from Ref. The conformal extrinsic curvature is 

decomposed into its trace K and tracelcss part Kij . The sign of the extrinsic curvature is that of Wald, with K > Q 
on a CMC hypersurface extending to 

The initial value equations on a CMC hypersurface constrain Kij through the momentum constraint. The well- 
known "zero-shear" condition necessary for to exist as a regular null hypersurface in the conformal manifold, 
expressed in our coordinates, says that 
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^B 
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Xb, 



k: 



Kl = 0. 



(8) 



The Hamiltonian constraint equation is an elliptic equation for fl, degenerate at where f2 = 0. The degeneracy 
allows no freedom in the first few terms of the expansion of the solution in powers of z. The asymptotic solution from 
Ref. H is 



n^^z 
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1 - ^^0^ 
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g ( ^0 + 7^0 - ^1 - X^^'XAB ) z 



C3 + ^Qlog^')z^ + 0(z^ 



(9) 



where C3 is a function of angle only known from the global (numerical) solution of the elliptic equation and 



-,AB 



X 



AB + Xb [ ^A - 1~^0Xa 



(10) 



The logarithmic term is present whenever outgoing radiation is present at It is a consequence of the CMC 

hypersurface condition. The coefficient Q also appears in the asymptotic expansion of the conformal extrinsic curvature 
derived from the momentum constraint equation: 



XbXaz 



di-Q log ■ 



Kz 



z^ + 0{z^). 



(11) 



If, in addition to the zero-shear condition, we require that the Wcyl tensor vanish at J^"*", the expansion of is 



-Xb+ 2^oXb2 + 0(z2) 



(12) 



and the expansion of K^ is 



^l = XATi3^ + ei^' + 0(z3). 



(13) 



The coefficients di and ei are functions of angle which are known only from a global solution of the momentum 
constraint equation. The vanishing of the Weyl tensor at ^+ , called the Penrose regularity condition in Ref. Q , is 
a necessary condition for the absence of polyhomogcneous terms in the Bondi-Sachs asymptotic expansion based on 
special null coordinates (see Ref. [isj). 

The time evolution of the conformal factor is given by 



dtil = P'^dk^ - i(5 ( K - OA' 



(14) 
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from which we see that K controls how the conformal gauge evolves. We want the conformal gauge to be consistent 
with being a two-sphere at the coordinate position z = at all times. Keeping ^+ at z = requires the boundary 
condition = ci = ao on the shift, so 9(51 = at z = 0. The only way to make sure that stays independent 
of angle is to require that dt£,o = 0. We also require fi^ = 0, in order that the angular coordinates be propagated 
along the null generators of J^^. These conditions greatly simplify the identification of the Bondi-Sachs mass aspect 
from the asymptotic behavior, since they make our coordinates close to inertial Bondi coordinates, for which the 
Bondi-Sachs mass aspect can be read off directly from the asymptotic expansion of the Bondi-Sachs metric. 

The evolution equation for ^, obtained from the evolution equation for the determinant of the angular part of the 
conformal metric in terms of the conformal extrinsic curvature, is 

a* (logo = P'dk{\ogO - + a Qlf - , (15) 

where f^^^c the divergence of the angular part of the shift considered as a vector field on the unit sphere. At ^+ , 
this becomes 

at(iogO = 5oQ^-ix). (16) 

Therefore, we impose the boundary condition 

K^l~Ko. (17) 

We make no other restrictions on the conformal gauge, i.e., on the details of the equation for K. As long as the 
equation is elliptic, as in Refs. Q and the boundary condition is easily enforced. If the conformal gauge is evolved 
as part of a purely hyperbolic system, in which there are no boundary conditions at things may not be so simple. 



III. THE BONDI-SACHS ENERGY AND MOMENTUM 



The Bondi-Sachs e nergy is the true physical energy defined at future null infinity in asymptotically flat spacetimes. 
Its original definition [14 . Il5| was based on a special class of retarded null coordinates (w, r, 6, with a metric of the 
form, in close to the notation of Ref. Q , 

ds^ = -Ve^^du^ - 2e^^dudr + r^fiAB {dx^^ - W^du) {dx^ - U^du) , 

where fiAB is a 2D metric asymptotically approaching the unit sphere metric Hab and whose determinant equals 
the determinant of hAB everywhere. We put bars on the angular coordinates to distinguish them from the angular 
coordinates in the CMC hypersurface. The expansion of Hab in powers of x = away from has the form 

hAB = hAB + XABX + ^X^'^ XCD flABX^ + 0{x^). (18) 

The absence of any traceless contribution to the 0{x'^) term in Eq. (jTS]) implies the vanishing of the Weyl tensor on 
y+, as Eq. (HD) did on CMC hypersurfaces. The remaining metric functions are asymptotically 

V = l-2M {u,e,^) x + 0{x^), (19) 
/3 = -^X^^XBca;2-HO(x3), (20) 

C/^ = -^X-'^TB^' + 0(x3)- (21) 

Like XAB, Xab {u, 0, (pj is a traceless symmetric tensor on the unit sphere. Its retarded time derivative along is 
the Bondi news. 

The Bondi-Sachs metric function Ad {ii^O^Cp) is what was originally identified as the Bondi mass aspect [l3l. Its 
average over the unit sphere at a fixed retarded time u is the Bondi-Sachs energy i^es a-nd the average over the unit 
sphere weighted by , where 

TV' = (sin 6* cos sin 6' sin (y9, cos 61), (22) 
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is the ith component of of the Bondi-Sachs momentum three-vector in the asymptotic Minkowski frame. However, 
there are other less coordinate-specific definitions of the mass aspect which give the same Bondi-Sachs energy and 
momentum when averaged. Some are based on the asymptotic behavior of Weyl tensor [111, [l^ HH US] and some are 
based on charge integrals derived from a Hamiltonian formalism (see, e.g., [l^jllSl)- What we will call the Bondi-Sachs 
mass aspect is not the metric function M, but rather the quantity which is monotonically decreasing on each of the 
null generators of [23I ]. 



Ma = M ~ 



AB 



lAB- 



(23) 



The monopole and dipole angular moments of A^a and M are identical for arbitrary . 

The task of actually calculating Ma from data on a CMC hypersurface can be accomplished in various ways. The 
charge integrals of Ref. [l7| require comparing the actual spatial metric and extrinsic curvature on the hypersurface 
with a background metric and extrinsic curvature. The results are sensitive to the choice of background quantities. 
The only fail-safe procedure mentioned in Ref. is to find the coordinate transformation from the CMC coordinates 
to Bondi coordinates, where the appropriate background spacetime metric is known, and then perform the inverse of 
this coordinate transformation on the Bondi-Sachs background metric. However, once one has transformed to Bondi 
coordinates, it is much simpler just to read off the Bondi-Sachs metric function M from the transformed metric. It is 
the latter procedure that we will follow. While the coordinate transformation involves the lapse and shift for the CMC 
coordinates, the final result only depends on the spatial metric and extrinsic curvature of a single CMC hypersurface. 

The lapse and shift which preserve the CMC slicing condition and the Gaussian normal coordinates based on 
were derived in Appendix A of Ref. @ . The singular terms in the elliptic equation for the conformal lapse allow no 
freedom in the leading terms of the expansion away from with the result 



-KqZ - 



0\ 



(24) 



The coordinate components of the shift are, using the boundary condition on K of Eq. (|17p . 



(25) 



and 



3^ = 5o?o' 



(26) 



Coordinate transformation to Bondi-Sachs coordinates 



The coordinate transformation to Bondi-Sachs coordinates is obtained from the transformation equations for the 
inverse metric. We replace r by its inverse x, so the Bondi-Sachs coordinates are (u, x, i'^) . The inverse Bondi-Sachs 
metric components, denoted by bars, are then 



0, 



x^Ve 



-2/3 -xA 



■2/3 



7,AB 



9 



(27) 



The metric transformation equations are in part partial differential equations for the Bondi-Sachs coordinates as 



functions of the CMC coordinates (t. 



and in part equations for the Bondi-Sachs metric functions. We solve 



these equations term by term in expansions in powers of z away from ^+ in the CMC hypersurface. From 



Tt) ^'Ttl)2^^[^z) ( 



du d 



& 
dz 



52 



{pr]+o{z-) 



= 0, 



we get 



^ = (« + r)^ + o(z^). 



The solution satisfying the boundary condition u ^ t in 

z 

u = t-\ 

2&0 



1 + -koz + 0{z'^) 



(28) 
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From 



uA 

9 = — 



oz oz ox^ 



+ 0(z3) 



= 



and 



du 



dx^ 8ao dx'-^ 



with the boundary condition x 



1 



-,AB 



i^-" ) + 0{z% 



The Bondi-Sachs x-coordinate is obtained from the transformation of the angular part of the inverse metric, 

,dx^ dx^ 



21.AB 



x'h 



\CD\ 



dx^ dx^ 



Equating the determinants of the two sides gives 



1 
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where A is the Laplacian operator on the unit sphere. Substituting back into Eq. pop gives 



1 X^^ 



-,AB 



X 



(29) 



(30) 



(31) 



(32) 



The second-order contribution to h^^ is consistent with the second-order contribution to h^^ once the correction is 
made for the time dependence of (Eq. ( A34) of Ref. Q ) , taking into account the difference between t and u at 
the same physical point in the interior of the liypersurface (Eq. (pS)) ). 

The Bondi-Sachs metric functions /3 and U"^ obtained from the coordinate transformation are consistent with the 
Bondi-Sachs asymptotic solutions from the Einstein equations, taking into account Eq. (|32[) . To obtain the metric 
function V , from which the mass aspect is derived, we use Eq. dn]) and the fact that 5 - = 0{z'^) to get 



x'^Ve- 



From Eqs. p4l [T5|) . the time derivative of Vl^ is 



(33) 



(34) 



Note that K cancels out, so there is no direct dependence on the conformal gauge. Inserting Eq. (|34|) into Eq. ([33 
and rearranging terms gives 



a 



(35) 



From the expansions of and ^ in Eqs. (jH |4]), 



= z 



1 1 , 

1 + -KQZ + g ( Co 



1 



Kz 1 



1 



c3+3giog^ + -eo% + ^^o^ 48 



6 12 

Combine Eqs. (I35|), ([Ml), Q, ([iT]), ([Ml), as weU as the first-order expansions of il, a, and /3^, and simplify to get 



(36) 



V = 1 



12 



nc + o{z^). 

(37) 
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The Bondi-Sachs mass aspect Ma as defined in Eq. ([23]) is 



Ma 



1 



1 



1 



-di + iQ ~ ^Co^Ko 



1 



2 8 6' 48 
with Q defined by Eq. pH)) . which is equivalent to 



(ki +xgxg) + 



1 



(38) 



(39) 



The ahernative approach to calculating the Bondi-Sachs mass aspect based on the Weyl scalar in the Newman- 
Pcnrose formalism is also straightforward to implement. The real part of ^2 is the component = s^s^Eki, where 

is the unit outward normal to the 2-surface, of the traceless electric part of the conformal Weyl tensor. Its 
normal derivative at ^+ is related to the mass aspect by 



Ma 



1 



t 2 "CD 

to X \CD 



M 



d.K.+X%dtxE 



(40) 



This simple form is valid if the Weyl tensor vanishes on .y^ and the gauge evolution boundary conditions preserve 
the explicit spherical geometry of J^^, with angular coordinates propagated along the null generators, such that the 
curves on ^+ with tangent vector dt are null geodesies. A 3D expression for Eij, from Appendix A of Ref. is 



TF 



(41) 



B. Conformally flat initial data 



The general expression derived in the previous section simplifies drastically for conformally flat initial data. By 
conformal flatness we mean that the conformal spatial metric is Euclidean, gij = Sij, but we still assume a CMC 
hypersurface. Future null inflnity is a coordinate sphere, 

R = V(a;i)2 + (a;2)2 + (a;3)2 = R^. 

The Gaussian normal coordinates based on are z = i?+ — R and the usual polar angles x"^ of spherical coordinates. 
The extrinsic curvature of constant- i? two-surfaces is isotropic, with k = 2/ R = 2/(i?_|- — z), and 

2 2 
^^0 = 7^ = 2^0, ^1 = 2^0 7 '^2 = 2^0 • 
R+ 

No gravitational waves are present at since xab = 'i'AB = 0. 

For the conformal extrinsic curvature on the CMC hypersurface, we adopt the generalized Bowen-York solutions of 
the vacuum momentum constraint equations from Ref. [l|: 

~ ~ C 3 

Q,~'^Kij = Aij [Snitij - 6ij] - [sikeS'^n^nj + SjuS'^n^ni] 

3 
3 

+ 7r75r + + '3''"'= (% ~ Sriirij)] , (42) 

with i?D |x - D|, n' = (x' - D') /i?D. Since the momentum constraint equations are linear, these solutions can be 
superimposed to generate initial data for multiple individual black holes or a single distorted black hole, depending 
on whether, after solving for the conformal factor, there are disjoint apparent horizons or a single all-encompassing 
apparent horizon. We follow Ref. [l| in only considering cases for which each solution represents a separate black 
hole centered at the coordinate position x^ = , and impose an inner excision boundary condition for each black 
hole so that the coordinate sphere i?D = ^^ms is a minimal 2-surface. The minimal surface condition, that the normal 
derivative of f2 vanish, is an inner boundary condition on the elliptic equation for the conformal factor. The choice of 
i?nis is an additional input parameter for each black hole. In the limit R-a^a — ^ 0, the Einstein-Rosen bridge associated 
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with the minimal surface becomes an infinitely long "trumpet" configuration. The minimal surface should be inside 
the apparent horizon, which means Rj^s should not be too large and also requires C > 0. 

For black holes on a conventional asymptotically flat (e.g. maximal) hypersurface, with the boundary condition 
— >■ 1 as i? oo, the sum of the boost 3- vectors P is the ADM linear momentum of the system. The physical 
interpretation of the boost is not so straightforward on CMC hypersurfaces. The vector Q represents the boost of 
the black hole as viewed from the other side of its Einstein-Rosen bridge. For a single centered black hole, it can be 
related to P through a condition of inversion symmetry about the minimal surface, but initial data without inversion 
symmetry are also valid. In our numerical results we will take Q = 0. Note that Aij is finite and non-zero at future 
null infinity, R = R^, which means that the traceless part of the physical extrinsic curvature goes to zero as fl^ for 
all contributions. This is quite unlike the situation on maximal hypersurfaces. 

Since conformal flatness implies that the conformal Ricci tensor vanishes, the equation for simplifies to 



nAn = 2 (^^'^^ —^^j^'' ' (43) 

with A the flat space Laplacian operator. The equation is elliptic almost everywhere, but is degenerate at the outer 
boundary R = where the boundary condition $7 = is imposed. Solutions are characterized by the values of K 
and i?+ and, for each black hole, the Bowcn-York parameters and i?ms- 

A natural dimensionless parameter for expanding the solution away from is z = ^o^- The asymptotic solution 
of Eq. (|^ then reduces to 



il = z 



1-^Z + CgZ^ (1 + Z) + O(Z^) 



(44) 



Here C3 = R\c3 is the rescaled dimensionless version of the locally undetermined coefficient C3 of Eq. The angular 
derivatives in Afl first affect the solution in 0{z^) and the conformal extrinsic curvature first contributes in 0{z^). 
With conformal flatness plus the regularity conditions, the 3D conformal extrinsic curvature has 

kl = n^N'AijN^ = diz^, 
and we define the rescaled dimensionless form of di as 

Ji = Rldi = 2Arr = (^] RlN'Ajm. (45) 
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Eq. ([55)1 for the Bondi-Sachs mass aspect simplifies to 



Ma = -^ (^4c3 + 1^1^ =-J^ (4c3 + Arb) . (46) 

Note that for spherical symmetry in Eq. (HU, all the Bowen-York parameters except C are zero and Arr — [K/if' C = 
C. See Ref. [l| for an extensive discussion of CMC hypersurfaces in the Schwarzschild geometry and the significance 
of the parameter C . The two contributions to the Bondi-Sachs mass aspect arc labeled by (Ma)^ for the contribution 
from the conformal factor and by (A/a)j^ for the contribution from the conformal extrinsic curvature. 

To sum up, the Bondi-Sachs energy and 3-monicntum associated with a CMC hypersurface are given by the integrals 

on 



and 



Sbs = — i MA\/det hAB dx^dx^ = — <t> MAsinedOdcj) (47) 
47r / * 47r / 

P^g ^ ^ <j> MaN"" sin eded<j). (48) 

An expression for the angular momentum as an integral over ^+ in the context of Bondi coordinates is given in Sec. 6 
of Ref. • In the absence of radiation, and specifically in the conformally flat case with Euclidean coordinates on a 
CMC hypersurface, this becomes 

J, = j) smOdedcjyeijkN' {r^A^''^ N*^. (49) 
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Analytic expressions for the (Ma)j^ contributions to the Bondi-Sachs energy and hnear momentum in terms of the 
Bowen-York parameters of Eq. (|42p are obtained in Appendix \^ and the angular momentum integral is evaluated in 
Appendix [B] The result for the angular momentum 3- vector J is deceptively Newtonian, 

J = S + DxP, (50) 

deceptive because D is a coordinate displacement, not a physical distance, and the boost P is not a physical linear 
momentum. 



IV. NUMERICAL EVALUATION OF THE BONDI-SACHS ENERGY AND MOMENTUM 

For conformally flat initial data, with an analytic solution of the conformal momentum constraint equation such as 
that of Eq. (j42[l . the only significant numerical task is solving the Hamiltonian constraint equation for the conformal 
factor ri. Our results are obtained with the pseudospectral elliptic solver that is part of the Caltcch-Cornell-CITA 
Spectral Einstein Code (SpEC), described in detail in Refs. [Tllll2[. The overall domain is divided into a number of 
sub-domains, which are spherical near the boundary at R ~ The output for the conformal factor in the outermost 
spherical sub-domain, from an inner radius i?i to is in the form of coefficients of a spherical harmonic expansion 
^ = J2e m ^imXl^ at each of a number of radial collocation points. 

To obtain the coefficients of the spherical harmonic expansion of C3, we fit the asymptotic form of Eq. to the 
numerical results for the conformal factor fi. It is counterproductive to use collocation points too close to R+ where 
the small contribution of the C3 terms to may not be large compared to numerical error, so we exclude collocation 
points with z < Z2 and fit 

to the form 

using a standard Grace^ curve-fitting routine. In practice, for single black holes, we find that taking — 0.95 

and Z2 > 0.02 gives values of (c3)qq stable to about 1 part in 10^, as long as the black hole centers are within 0.12 i?+ 
or so, and the boost parameters are not too large. We also need the coefficients of the spherical harmonic expansion 
of Afm, which are obtained analytically in Appendix |X] for I ~ and t ~ 1. Then the Bondi-Sachs energy is 

£;bs = -| [4 (c-3)oo + {Ab.b)J \[^. (51) 

The Cartesian components of the radial unit vector are composed only oi t ^ \ spherical harmonics, and since 
C3 is real, the standard Condon-Shortley phase convention gives {c^)i_i = ^{^3)11- Then from Eq. (|48|) for the 
Cartesian components of the Bondi-Sachs momentum, 

Pbs = +|R44(c3)n + (^i?i^)n] 

^B^s = -|lnr [4(c-3)n + {Ani,) (52) 

^BS = -|[4(c3)io + Mj/^- 

Analytic expressions for the contributions to the Bondi-Sachs energy and 3-momentum from the conformal extrinsic 
curvature, the Aj^ji terms in Eqs. (|5ip and (j52p . and their dependence on the Bowen-York parameters, are derived 
in Appendix [X] The results are given in Eqs. (jA2[) and (|A4p . The total Bondi-Sachs energy and 3-momentum form a 
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4-vector in the asymptotic Minkowski spacetime, and the Bondi-Sachs mass of the system is defined to be the special 
relativistie magnitude of this 4- vector, 

Mbs = V^^Is-Pls- (53) 

The angular momentum of the system has no contribution from the conformal factor and is given analytically by 
Eq. dnOl) (see Appendix H]). 

There are two arbitrary scale factors in our problem, the choice of the unit of mass-length-time and in the definition 
of the conformal factor. Only combinations of variables invariant under both a change of units and a uniform conformal 
rescaling are physically meaningful. Combinations of the input parameters with this property are: 

C = {K/ifC, = {K^R+/l^)P\ = {Kyi8R+)Q\ = S'/C, R^s = Rms/R+, = D'/R+. (54) 

For sufficiently small i?ms; so there is a trumpet-like configuration inside the apparent horizon. ~ {K/3)M for a 
Schwarzschild-like black hole of mass AI (see Fig. 4 in Ref. [H). Making C < 1 ensures that the warping effects of 
the mean extrinsic curvature are small in the vicinity of the apparent horizon of a centered black hole. As argued in 
Ref. [l[, when this is the case, the physical momentum should be roughly fimax times the boost. From the asymptotic 
solution for the conformal factor, fimax ~ \d^/9R\fi^ (R+/2) = KR+/6, so that we expect KR^/6 times the Bowen- 

York boost (or 3/K times P*) to be a rough estimate of the physical Bondi-Sachs 3-momentum of the black hole. 
Inversion symmetry about the apparent horizon, as sometimes assumed, implies = ±i?^gP% so at the terms 
in the conformal extrinsic curvature are of order R^^ times the terms and typically make a negligible contribution 
to the Bondi-Sachs momentum. We will simply take = in our numerical calculations. Well away from any black 
hole apparent horizons and not too close to P_|_ (R < P_|_/3), the conformal factor plateaus near its maximum value, 
so in this regime physical distances are roughly coordinate distances divided by KRj^^/Q. 

A key output parameter is the dimensionless rescaled "irreducible mass" Mi„ = {K/i)M;,, = (A73)^^/167r, where 
A is the area of the black hole apparent horizon. The irreducible mass is used to scale the results for the Bondi-Sachs 
energy, momentum, and mass. In Ref. the irreducible mass is used as a surrogate for the physical mass in the 
discussion of results, so giving our results in units of the irreducible mass facilitates comparison with that paper. 

For single black holes, we vary the Bowen-York parameters starting from a representative spherically symmetric 
model that has C = 0.0011207, R^as = 0.00127, and = 0.028339. For these parameters, the effects of the non-zero 
K start becoming important at about 10 times the radius of the apparent horizon, which means initial data near 
the black hole are similar to initial data on a conventional maximal hypersurface. We also present results for two 
examples of binary black hole initial data. 

A. Single spinning black hole 

Here we see how the Bondi-Sachs mass for a single centered black hole with no Bowen-York boost varies with 

Bowen-York spin S = (0, 0, S). As the dimensionless magnitude of the spin S increases, we keep CcS = C\J 1 + 2 (S*)^ 

rather than C constant, because the square root of CeS tracks Mi„ much better than as S becomes larger than 1. 
There is no linear momentum, so the Bondi-Sachs mass equals the Bondi-Sachs energy, and by Eq. ([50| the angular 
momentum equals the spin. The ratio MBs/-^irr — 1 is plotted in Fig. [T]as a function of the spin extremality parameter 
C = 5/(2Mj^,j.), for < C < 0.78. The sequence terminates at C = 0.78 in the hmit S oo for eonformally flat initial 
data (see Ref. [l|). We also plot the corresponding quantity for the Christodoulou mass Mq, which is defined as 
the mass of the Kerr black hole with the same Mjii. and angular momentum. The Bondi-Sachs mass for our initial 
data is a bit larger than the mass of the Kerr black hole for the same Q, but the maximum Afc/Min. — 1 is significantly 
larger, 0.414, because for a Kerr black hole the maximum ^ and the maximum J /Mq arc both 1. 

B. Single displaced black hole with zero spin and boost 

Now we displace a single, non-spinning, unboosted black hole a coordinate distance D (in the x-direction) from 
the center of the coordinate sphere. The other Bowen-York parameters are the same as for our representative 
spherically symmetric black hole. Surprisingly, rather than being zero, the velocity of the black hole in the asymptotic 
Minkowski frame (given by the ratio of the Bondi-Sachs momentum to the Bondi-Sachs energy) is proportional to 
and in the same direction as the displacement from the origin, as shown in Fig. [21 The dominant contributions to the 
Bondi-Sachs energy and momentum are from the conformal factor, and have opposite signs from the contributions 
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FIG. 1: Single spinning unboosted centered black hole: The ratio of the Bondi-Sachs mass (Mbs) to the irreducible mass (Mirr) 
is compared to the ratio of the Christodoulou mass (Mc) to the irreducible mass, as the spin extremality parameter ( is varied 
from to 0.78. 




FIG. 2: Single black hole displaced in the i-direction with no boost or spin. The velocity of the black hole in the asymptotic 
Minkowski frame, \Pgg\/EBS, is plotted against the relative displacement = D^/R+. The remaining input parameters {C, 
Rms) have their usual values. 



associated directly vifith the extrinsic curvature. The contributions from the extrinsic curvature, as calculated using 
Eqs. ^ and jSl), are simply (i;Bs)K = -KC/3 and (Pbs)k = -(D/i?+)ifC/3. With D = 0.1, Mi„. = 0.028340, 
and C — 0.0011207, the fractional extrinsic curvature contributions to the Bondi-Sachs energy and momentum are, 
respectively, (£;bs)k/£^bs = -0.039 and (Pbs)k/Pbs = -0.019. 

As an indication of whether the displaced black hole is just a moving Schwarzschild black hole or whether it has 
additional distortions, we compare the Bondi-Sachs mass as given by Eq. (|53|) with the irreducible mass derived from 
the area of the apparent horizon. We find that these are the same to the numerical accuracy of our calculations 
(about 5 significant figures), which suggests that the displaced black hole is in fact just a Schwarzschild black hole in 
an asymptotic Lorentz frame in which the black hole is not at rest. We have not tried to prove this more decisively, 
by finding the coordinate transformation to the conventional form of the Schwarzschild metric. 
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C. Single centered boosted black hole 

Now consider a Bowen-York boost as the only deviation from Scliwarzschild. On an asymptotically flat hypersurface, 
such as a maximal hypersurface, the Bowen-York boost is equal to the physical ADM momentum. In the limit of 
small K^C, there is a range of R over which is roughly constant near its maximum and Mj^r << R/^ << ^/K, 
a region of the CMC hypersurface approximating the asymptotic part of a maximal hypersurface. Ref. [l[ argued on 
this basis for using P riniax an estimate of the Bondi-Sachs linear momentum. The close agreement between and 
(i^/3)P]|g in Tabic U shows that over the entire range of boosts considered, {KR^/6)P^ is within about one percent 
of P^s- The ratio of the actual numerical result for r^max to the KR+/6 analytic estimate varies from 0.931 to 0.858 
over this range, so it is considerably better to use the analytic estimate for f2max, rather than its actual numerical 
value, in relating the Bowen-York boost to the physical momentum of a centered black hole. 

The fact that MesZ-^^irr becomes significantly larger than 1 as the black hole motion becomes relativistic shows 
that the boosted Bowen-York black holes are definitely not simply Schwarzschild black holes moving relative to the 
asymptotic Minkowski frame. There is additional energy associated with non-kinematic distortions of the geometry. 
From Eq. (|A2p . we see that for centered black holes, the direct contribution to the Bondi-Sachs energy from the 
conformal extrinsic curvature is independent of the boost, and for the initial data calculated here, (i?Bs)K/-^irr ~ 
—0.04, which is small compared with the positive contribution from conformal factor. In contrast, the the Bondi-Sachs 
momentum is strongly dominated by the contribution from the conformal extrinsic curvature for centered black holes, 
since by Eq. dH, (Pbs)k = (3/i^)P. 

The boosts are taken to be zero rather than the inversion-symmetric values often assumed in Bowen-York 
initial value calculations on maximal hypersurfaces (e.g., in Ref. [25|). Only one side of the Einstein- Rosen bridge 
is astrophysically relevant, since the other side cannot exert any causal influence on the astrophysical side. There is 
no very compelling reason, other than geometrical elegance, to assume any symmetry between the two sides of the 
bridge or that the bridges of multiple black holes connect to a common asymptotically flat space. The geometry of 
the physical space in the immediate vicinity of the black holes may depend to some extent on whether the are 
zero or inversion symmetric, with possible small effects on the results for the Bondi energy and momentum, but such 
effects should not be noticeable for the near trumpet considered here, and disappear completely in the limit of an 
extreme trumpet configuration {Rms 0). 
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2.32632 


3.39577 
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5.67844 
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1.0 


1.00261 


1.03387 


1.28624 


1.74330 


2.30348 


2.92016 



TABLE I: Centered black hole with varying Bowen-York boost . The C parameter is fixed at our adopted standard value, 
and Rms/R+ is adjusted to keep the irreducible mass roughly constant as P^ increases. The input and output quantities are 
tabulated in the dimensionless and conformal-scale invariant forms discussed earlier in Section HVl 

For a slightly different slant on our results, we plot in Fig. [3] (i) the fractional "excess mass" of the black hole 
-^BsM^^irr — 1, (ii) the velocity of the black hole = P^q/Ebs in the asymptotic Minkowski frame, and (iii) the 
similar excess ADM mass for a boosted Bowen-York black hole on a maximal hypersurface, all as functions of the 
ratio of the appropriate physical momentum (Bondi-Sachs or ADM) to the irreducible mass. Data for the last curve 
are taken from Table II of Ref. [1^, for C = and inversion-symmetric Q'. Note that as the boost and physical 
momentum increase, the gravitational mass Mbs (solid black curve) increases rapidly enough that the black hole's 
velocity = Pgg/E'BS (violet dash-dot curve) asymptotes to a value less than one (about 0.86). The results of 
Ref. [25| are quite similar, but with a slightly higher asymptotic velocity. The excess gravitational mass can be 
thought of as energy stored in distortions of the geometry relative to a boosted Schwarzschild black hole, energy 
which would be radiated into the black hole and to infinity during subsequent evolution. 

The black holes on CMC hypersurfaces have more excess mass than the corresponding black holes on maximal 
hypersurfaces, modestly so as long as the motion of the black holes is non-relativistic, but more substantially as the 
motion becomes relativistic. The Cook and York initial data has C = 0, which in the case of spherical symmetry 
means that their maximal slice passes through the intersection of the future and past Schwarzschild event horizons, 
as opposed to the more trumpet-like configurations chosen here and emphasized in Ref. [l[ . How much effect this has 
on the difference in behavior needs further calculations to sort out. 
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FIG. 3: Centered black hole with Bowen-York boosts. Using data from Table U above and Table II of Cook and York [25l ]. 
we plot (i) the fractional excess of the Bondi-Sachs mass over the irreducible mass (solid black line), (ii) the corresponding 
excess ADM mass of Bowen-York initial data on maximal hypersurfaces (dashed red curve) , and (iii) the velocity of the black 
hole — P^s/Ebs (dash-dot violet line) in our initial data. All quantities are shown as a function of the ratio of physical 
momentum, Pbs or Padm as appropriate, to the irreducible mass. The curves are interpolations between results for discrete 
models. 

D. Displaced and boosted black holes 

In this section, we lay the groundwork for building a binary black hole initial data set on CMC hypersurfaces. 
For binary black holes in circular orbit, we want each black hole considered by itself to have zero radial momentum 
and transverse momentum corresponding to that of a mass in a circular orbit, which can be estimated crudely from 
the Newtonian equations of motion. Since we have seen that a displaced unboosted black hole on a CMC slice has 
a non-zero radial momentum, this requires finding the radial boost which makes the net Bondi-Sachs momentum 
of a displaced single black hole zero, and then the transverse boost which generates the desired transverse physical 
momentum. Once this is done for each black hole of the binary considered separately, the Bowen-York data for the 
black holes can be superimposed and the Bondi-Sachs energy and momentum of the full system calculated. As long 
as the two black holes are well-separated, we should find that the Bondi-Sachs energy of the combined system is less 
than the sum of the energies of the black holes considered separately by something close to the Newtonian binding 
energy of the system, provided that the total Bondi-Sachs momentum of the system is close to zero. 

The first step is to construct a single black hole which is displaced from the center of the coordinate grid and at 
rest. Fig. |4] shows that this result can be achieved by giving the black hole an appropriate boost antiparallel to its 
displacement. For the example shown in Fig. 21 the hole is displaced in the positive x-direction, with = 0.1. A 
rescaled boost = —0.006 very nearly cancels the radial momentum associated with the displacement, yielding a 
net P|g/Mirr = -9.6 x 10"^. 

Next, we build an initial data set with two non-spinning displaced black holes, each of which considered separately 
has zero Bondi-Sachs momentum, and calculate the gravitational binding energy; i.e., the decrease in the Bondi-Sachs 
energy of the combined system compared with the sum of Bondi-Sachs energies of each hole separately. The input 
parameters are chosen so that, as estimated by the procedures outlined in Sec. HID of |l| for Schwarzschild black 
holes, the ratio of minimal surface radius to horizon radius for each hole is about 0.8, the ratio of the black hole masses 
M^/M^ 2/1, and the total rescaled irreducible mass of the system is about 0.7 that of our standard single black 
hole. We locate holes A and B along the x-axis at Z)^ = 0.04 and = —0.08, so the center of mass of the system 
is close to the origin. The values for and are, respectively, 0.000228 and 0.0000568. The scaled minimal 
surface radii are = 0.000978 and = 0.000489. We apply Bowen-York boosts which for each black hole in 
isolation approximately cancel the Bondi-Sachs radial momentum due to that hole's displacement. The dimensionless 
Bondi-Sachs masses {Mbs = {K/3)Mbs) for each hole in isolation are M^g = 0.013104 and M^g = 0.006604. When 
the black holes are superimposed, we find that the combined Bondi-Sachs momentum is less than the sum of the 
already small values for the individual holes. The dimensionless Bondi-Sachs mass of the superimposed black holes 
is Mbs = 0.019368, which is 0.000340 less than the sum of the individual values. The Newtonian gravitational 
binding energy of the system is U = M^M^ /tab ~ ^ma^M^M^ /ID"^ — or, using the analytic estimate of 
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FIG. 4: Finding the Bowen-York boost that cancels the Bondi-Sachs momentum associated with the displacement — 0.1. 
The CMC parameters have their usual values, except that is varied with the boost to keep Mirr approximately constant. 



^^max, KU /3 ~ 0.5-R+ -^Bs -^Bs/l^^ ~ "^^1 = 0.000371, which is reasonably close to the rcscalcd general relativistic 
binding energy. 

The third step is to determine how Bowen-York boosts transverse to a black hole's displacement translate into the 
transverse motion of the black hole as revealed by its physical transverse momentum. We have calculated the Bondi- 
Sachs energy, linear momentum, and angular momentum of a displaced black hole with nearly zero radial momentum 
and various transverse boosts. The relationship between transverse boost and transverse physical momentum is very 
close to what we found between boost and momentum for a centered black hole in Sec. IIV CI An example of the 
comparison for a black hole given a moderately relativistic Bowen-York boost transverse to the displacement in the 
x-direction is given in Table HIl We conclude that it is acceptable to use the results for a centered boosted black hole 
in Sec. IIV Cl including the analytic formula, when estimating the transverse Bowen-York boosts appropriate for black 
holes in quasi-circular binary orbits. Of course, the displaced black holes also acquire orbital angular momentum 
according to Eq. ([50)) . 
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0.13889 
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0.13889 
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0.47686 


1.11117 


1.00365 


3.45956 



TABLE II: A centered black hole (first row) and a displaced black hole (second row) are both given the same mildly relativistic 
(transverse) rescaled Bowen-York boosts P^ — 0.013889. The displaced black hole is also given a longitudinal rescaled boost 
P^ = —0.00660 to approximately cancel the physical momentum associated with the displacement. The other input parameters 
for both holes have the standard values. The black holes have almost identical physical properties, except for the orbital angular 
momentum of the displaced hole. 

The close relationship between the transverse Bondi-Sachs momentum and the transverse rescaled Bowen-York 
boost can be understood as a consequence of the dominance of the direct contribution to the Bondi-Sachs momentum 
from the conformal extrinsic curvature, as calculated in Eq. (|A4|) . over the contribution from the conformal factor 
for these transverse components. As long as D is reasonably small, the dominant transverse term in Eq. (|A4[) is just 
(KR^/6)Ptra,ns- There is no simple relation between the longitudinal component of the Bondi-Sachs momentum and 
the longitudinal component of the Bowen-York boost because the conformal factor contribution to the longitudinal 
component is not small. 



E. Binary black hole initial data 



Here we build on the results of the previous subsections to construct initial data for two spinning unequal-mass 
black holes with initial momenta roughly appropriate for circular orbits. The input parameters are chosen based on a 
Newtonian analysis, in circumstances such that the post-Newtonian corrections are roughly at the 10% level. Ignoring 
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relativistic corrections is consistent with other sources of uncertainty, such as estimating parameters for each black 
hole in isolation, ignoring their tidal interactions. While similar initial data were presented in [l|, we are now on much 
firmer ground in considering the relationship between the Bowen-York parameters and the physical properties of each 
hole. 

We consider two examples. Both are constructed to have about a 2:1 mass ratio and roughly a 12:1 ratio of physical 
distance between the black holes to the total gravitational mass of the system. The minimal surface radii are chosen 
to give moderately trumpet-like configurations. The black holes are assigned close to maximal spins oriented in the 
orbital plane and perpendicular to the displacement between the holes. In one case (Case I), the ratio of coordinate 
separation of the holes to the coordinate radius of J"^ is 0.12, and in the other case (Case II), it is three times larger. 
In the latter case, there is the potential for greater computational efficiency, since the black holes can be resolved with 
a lower resolution grid. 

The Newtonian analysis of a circular orbit binary system with masses Ma and Mb separated by a distance tab 
gives for the equal and opposite momenta of the masses 



perpendicular to the displacements of the black holes (which are along the x-axis). Equating the transverse rescaled 
Bowen-York boost to KjZ times the transverse physical momentum, for A/a ~ 2Mb, we choose = 0.064(7^/3) (Ma -I- 
Mb), since tab is about 12 times the total mass. Estimating black hole masses from the Bowen-York parameters is a 
bit uncertain when the black holes are not Schwarzschild, but consistency with the total system mass can be verified 
retrospectively. 

The non-zero input parameters chosen for Case I are Ca = 2.4531 x 10-^ Cb = 5.1224 x 10-^ S\_= 6.5224, 5*1 = 



-7.8088, D\ = 0.04, = -0.0_8, i?4s = 9.78 x lO""*, B^^^^ = 4.8 9 x 10 '^, and boost components PI = -0.001074, 
P\ = 0.00116, Pjf = 0.001060, P| = -0.00116. As noted in Sec. HvXl we expect the rescaled black hole masses to 



of iS* for each hole mean that the black holes have close to the maximum possible spin extremality parameter C,. 
_In Case II, _the input parameters are Ca = 3.6231 x 10-^ Cb = 8.5343 x 10-^ S\ = 3.9347_, = -4.1761, 
U\ = 0.12, = -0.24, i?^^, = 2.934 x lO'^, B^^^ = 1.467 x lO'^, and boost components P| = -0.010026, 
PI = 0.00348, P| = 0.010617, P| = -0.003639. The Coff values are 0.04526 for hole A and 0.02261 for hole B, about 
three times the values in Case I, while the D values are also three times larger. This suggests the two cases are similar 
in the ratio of physical separation of the holes to the total mass of the system. The S values are smaller in Case II 
than in Case I; this is because the elliptic solver failed to converge if we tried to make the S values as large as in Case 
I. Another effect of the larger ratio of the size of the binary system to ZjK in Case II is the larger P^ values needed 
to cancel out the radial components of the individual black hole momenta. 

Our results for the irreducible masses and angular momenta of the individual holes and the physical energy, mo- 
mentum, mass, and angular momentum of the systems are given in Table IIIII Note that there are non-zero, though 
small, z-components for the system velocities even though none of the boosts has a non-zero z-component. This is 
the result of the spins of the black holes interacting with the geometry of the hypersurface, one indication of which 
is the D x S term in Eq. (|A4[) . Also, note that Kj?) = 0.02 gives a total system mass close to 1 in Case I, and this 
together with P+ = 100 makes coordinate distances close to physical distances. The same is true for Case II since K 
is increased and P+ is decreased by a factor of 3. The calculations for Case II took about 80% of the time on the same 
computer as the calculations for Case I, indicating some potential for realizing increased efficiency of computation by 
having the binary system occupy a larger fraction of the conformal domain. 

While the system momentum and velocity should ideally be zero, this is not precisely the case in our calculations, 
which made no attempt to fine-tune the input parameters. Still, the residual system velocities are quite small compared 
with the Newtonian estimates for the orbital velocities. The system mass is significantly greater than the sum of the 
irreducible masses in both Cases I and Case II, since the contributions of the spins to the black hole hole masses 
dominates over their gravitational binding energy. The slightly smaller spins in Case II make this effect somewhat 
smaller there. Of course, the system mass can never be less than the total irreducible mass. 




scale roughly as Ccff = C^J 1 + 2(3^ , which is 0.0150864 for hole A and 0.0075366 for hole B. The large values 



V. 



DISCUSSION 



The first part of this paper is a discussion of the asymptotic behavior of the conformal factor fi, the conformal spatial 
metric, and the conformal extrinsic curvature of CMC hypersurfaces in the neighborhood of .y^ as constrained by the 
initial value equations for astrophysically general asymptotically flat spacetimes, and a calculation of the Bondi- Sachs 
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Case I 




















hole A: 


1.06367 


n/a 


n/a 


n/a 


n/a 


n/a 





1.60000 


0.92800 


hole B: 


0.53449 


n/a 


n/a 


n/a 


n/a 


n/a 





-0.40000 


1.85600 


system: 


1.59816 


1.96685 


1.96685 


-1.04379 


0.25464 


0.05850 





1.20000 


2.78400 


Case II 




















hole A: 


3.16247 


n/a 


n/a 


n/a 


n/a 


n/a 





14.25600 


8.35200 


hole B: 


1.62032 


n/a 


n/a 


n/a 


n/a 


n/a 





-3.56400 


17.46720 


system: 


4.78280 


5.90011 


5.90003 


-5.25838 


-0.16546 


0.17739 





10.69200 


25.81920 



TABLE III: Physical quantities for two binary black hole initial data sets. (See text for input parameter values.) From left 
to right are the irreducible masses, the Bondi-Sachs energy, the Bondi-Sachs mass, the system velocity components, and the 
angular momentum components, all made dimensionless by appropriate powers of K/?>. There are no well-defined Bondi-Sachs 
energy and linear momentum for the individual black holes, but the angular momenta can be calculated using Eq. (|B7|) . 



mass aspect in terms of the coefficients in the expansion away from the ^+ 2-surface, which is the intersection of the 
CMC hypersurface with the ^+ null hypersurface. In this analysis, we adopt special spatial coordinates, Gaussian 
normal coordinates based on and use the conformal proper distance z from along the normal spatial geodesies 
as the expansion parameter. Angular coordinates are propagated along the normal spatial geodesies, initialized to be 
standard polar coordinates on which is assumed to be a true sphere with intrinsic curvature ^q. This is not a 
physical constraint for an asymptotically flat spacetime with a topologically spherical null infinityjiust a constraint on 
the conformal gauge. The asymptotic expansion is essentially identical to that presented in Ref. [9[, except we do not 
allow the polyhomogeneous terms in the expansion of the angular part of the conformal spatial metric through 0{z'^) 
that are a major focus of that paper. No such terms are present unless they are present at some arbitrarily early time, 
which means that they should be associated with incoming radiation at past null infinity. In typical astrophysical 
problems incoming radiation near ^+ should only be generated by backscatter or nonlinear interaction of outgoing 
radiation, and should not be present ah initio, except for what might be plausibly associated with backscatter of 
radiation emitted before the start of the numerical calculation. 

As noted in Rcf. Q, polyhomogeneous terms are nevertheless gcnerically present at 0{z'^) in the expansion of the 
conformal factor and at 0{z^) in the expansion of the conformal extrinsic curvature. The authors of Ref. Q considered 
the vanishing of the expression in Eq. ((39)) a condition for the absence of these polyhomogeneous terms, but we point 
out that the condition cannot be satisfied over any finite time interval in which outgoing gravitational radiation, as 
indicated by time variation of xab, is present at J^'^ (see also Ref. Q). The implied smoothness of only for the 
conformal factor and for the conformal extrinsic curvature at J^~^ is a gauge artifact of the CMC hypersurface 
condition. The offending terms are not present if the physical mean curvature K has precisely the right angular 
dependence at 0{z^) in a power series expansion away from This cannot be accomplished by any predetermined 
specification of K, but maintaining a high degree of smoothness at J^"*" is possible with fully hyperbolic evolution 
schemes, such as the regular conformal field equations of Friedrich [1^, H^l, as well as with the Bondi-Sachs gauge. 

For the assumed special geometric properties of to be preserved during evolution, a constraint must be placed 
on the evolution of the conformal factor; specifically, that the trace of the conformal extrinsic curvature K satisfy the 
boundary condition of Eq. (jl7|) . This says that the 3D mean conformal extrinsic curvature of the CMC hypersurface 
at which governs the evolution of the conformal factor, should equal the 2D mean extrinsic curvature of the 
^+ 2-surfacc when both are rcnormalized according to their respective number of dimensions. All this also assumes 
preservation of the coordinate location of ^+ as a coordinate sphere centered at the origin and propagation of the 
angular coordinates on along its null generators, which imposes boundary conditions on the shift vector. These 
boundary conditions can be implemented in a mixed hyperbolic-elliptic evolution scheme like those of Ref. Q or 
Ref. Also, we impose the usual regularity conditions at J^~^ , the "zero-shear" condition and the vanishing of the 
Weyl tensor, which as shown for the evolution schemes of 0, [1], are preserved by the evolution equations. 

We have applied this asymptotic analysis to derive the expression for the Bondi-Sachs mass aspect given in Eq. ([55]) . 
from which the Bondi-Sachs energy and momentum can be calculated as appropriate angular averages over This 
is a 3D expression, involving only quantities evaluated on a single CMC hypersurface. The 3D expression derived in 



Section 5.3 of Ref. jI7l | requires knowledge of a "background" metric and extrinsic curvature in addition to the actual 



metric and extrinsic curvature of the CMC hypersurface. It is non-trivial to determine the appropriate background 
quantities. The terms appearing in Eq. arc simply related to the asymptotic behavior of the conformal factor. 
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conformal metric, and conformal extrinsic curvature in Gaussian normal coordinates, and can be determined by solving 
the geodesic equation for the geodesies normal to ^+ in the CMC hypersurface in any computational spatial gauge 
satisfying our boundary conditions. The alternative expression of Eq. (|40)) based on the Weyl tensor is completely 
equivalent, as can be seen by evaluating the Weyl tensor in the form of Eq. (|4ip in Gaussian normal coordinates. 

With conformally flat initial data, all of this simplifies drastically. With J'^ a coordinate sphere of constant radius 
in a flat conformal geometry, the normal geodesies are radial, and the 2D extrinsic curvature of surfaces of constant 
coordinate radius R is just times the unit sphere metric. The momentum constraint equation admits solutions 
identical in form to the Bowen-York solutions. We have shown in the Appendices how to evaluate analytically the 
angular integrals of the term coming from the conformal extrinsic curvature in the mass aspect in terms of the 
Bowen-York parameters. The only numerical calculation necessary is solving the Hamiltonian constraint equation for 
the conformal factor and fltting to a power series expansion in the coordinate distance from taking advantage 
of the known analytic form through order . We have computed how the Bondi-Sachs energy- momentum and the 
angular momentum of the system depend on the Bowen-York parameters for a sampling of single and double black 
hole systems, similar to those considered in Ref. [l|. 

For a single black hole centered at the origin, with varying Bowen-York spin (which in this case equals the system 
angular momentum), the results are displayed in Fig. [T] They show that the "rotational energy" of the black hole, 
the excess of the Bondi-Sachs mass over the irreducible mass, has a dependence on the spin extremality parameter 
C similar to that of Kerr black holes. The Bondi-Sachs mass is a bit larger in our case, as one would expect. The 
main difference is that our sequence cannot be extended to as large a value of the extremality parameter as the Kerr 
sequence. 

A single non-spinning, unboosted black hole displaced from the origin, i.e., displaced relative to the center of the 
sphere representing future null infinity, is more interesting. As shown in Fig. [21 the displaced black hole acquires 
a velocity in the asymptotic Minkowski frame in the direction of its displacement close to twice the ratio of the 
displacement to Consequently, when constructing initial data for black hole binaries on CMC hypersurfaces, 

compensating Bowen-York boosts directed inward should be introduced for each displaced black hole, so that its net 
radial velocity is zero. 

Results for a single centered black hole given various, mostly rather large, boosts are presented in Sec. IIVBI Here 
we emphasize how close the numerical result for the Bondi-Sachs momentum Pbs is to KR^ /6 times the Bowen- 
York boost vector, even when the black hole is moving at relativistic speeds. Ref. (l|, presented a crude hand- waving 
argument that the physical momentum should be roughly fimax (calculated numerically) times the boost. The expected 
accuracy of this estimate was of order KMi„, or about 10% for our choices of Bowen-York parameters. Comparison 
of the first two rows of Table U shows that times the boost is within about 1% of Pbs, while differs 

from flniax by about 10%. The physical properties of these boosted centered black holes are explored further in Fig. [3] 
and compared with corresponding results for boosted black holes on maximal hypersurfaces. In both cases, the ratio 
of physical mass to irreducible mass increases rather rapidly with increasing boost, and the velocity of the black hole 
plateaus at a value somewhat less than the speed of light for large boosts. Bowen-York initial data cannot produce 
black holes moving at speeds arbitrarily close to the speed of light. 

A combination of boosts and displacements is necessary to construct initial data for binary black holes. In Sec- 
tion IIV Dl we first determine the radial boost opposite to the direction of the displacement necessary to cancel the 
radial momentum induced by the displacement. Orbital motion of a black hole in a binary system requires a transverse 
boost, and we verified in Table HIl that the relationship between the transverse boost and transverse momentum for 
a displaced black hole is very nearly the same as the relation between boost and momentum for a centered boosted 
black hole, as calculated in Section HVBl We also checked that the gravitational binding energy of two black holes at 
rest separated by about 10 times the total gravitational mass of the system agrees reasonably well with the Newtonian 
formula. 

Finally, we make a crude attempt to construct initial data for two black holes in circular orbits about the center 
of mass. Again, the black holes, with about a 2:1 mass ratio, are separated by a proper distance about 10 times the 
gravitational mass of the system. Boosts are assigned so that, considered individually, each black hole has zero radial 
momentum and transverse momentum roughly as expected from a Newtonian analysis of the orbit. The black holes 
are also given close to maximal spins oriented in the plane of the orbit. The total energy, momentum, and angular 
momentum of the system are calculated and displayed in Table IIIII The numerical results are consistent with the 
system linear momentum being close to zero, in that the system -Pbs/^bs values are rather small compared to the 
P'/Mirr for each black hole. The component of the residual momentum along the line of the black hole displacements 
(the x-component) is considerably larger than the other components, suggesting that we should have taken account of 
the black hole spins in estimating the radial boosts. The system angular momentum, including both spin and orbital 
angular momentum, is roughly consistent with Newtonian expectations. The two cases we present are designed to 
have similar ratios of physical separation to black hole masses, but values of {K/3)Mbs differing by about a factor 
of three. The ratio of the coordinate separation of the holes to R+ is three times larger in Case II than in Case I, 
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implying correspondingly larger effects of the warping of the CMC hypcrsurface at the scale of the binary system. 
Nonetheless, Case II achieves close to the same accuracy as Case I with half the number of collocation points in the 
outer spherical computational domain surrounding the two holes, 40 versus 80. 

The exploratory analyses of conformally flat initial data on CMC hypersurfaces carried out in Ref. [ij and this paper 
need to be extended in various ways to approach the current sophistication of initial data calculations on maximal 
hypersurfaces. In particular, it would be interesting to see how well the numerical techniques used here hold up when 
considering non-conformally-flat data, particularly when radiation is present at How well does the SpEC elliptic 
solver handle the more complicated asymptotic behavior, including the polyhomogeneous terms in the conformal 
factor and extrinsic curvature? This is important for implementing mixed hyperbolic-elliptic evolution systems. Also, 
quasi-equilibrium methods [28| - |30| based on the conformal-thin-sandwich approach to the initial data problem [3l|, |32| 
have come to the fore in recent years as a way of obtaining quieter starts to binary black hole calculations, with less 
junk radiation and more circular initial orbits. It is clearly worthwhile to attempt to adapt these approaches to initial 
data on CMC hypersurfaces. Can initial data for more Kerr-like black holes, as in Ref. [Sj], be constructed on CMC 
hypersurfaces? We hope to explore such issues in future papers. 
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Appendix A: Analytic calculation of the partial Bondi-Sachs four-momentum from the conformal extrinsic 

curvature 



The Bondi-Sachs mass aspect for the conformally flat initial data is given in Eq. (|46)) as the sum of two terms, 
one a coefficient in the expansion of the conformal factor away from ,^+, and the other a projection Ann of the 
conformal extrinsic curvature defined by Eq. (j45|) . The first term requires solving a non- linear elliptic equation, 
which can only be done numerically. However, in the case of Bowen-York initial data on a conformally flat spatial 
geometry, the conformal extrinsic curvature Aij is as given explicitly in Eq. (|42p . The angular integrals on ^+ for 
the Ajiji contribution to the Bondi-Sachs mass and linear momentum will be obtained in analytic form with the help 
of Mathematica. 

For a given black hole, we can rotate the Cartesian coordinates in the conformal flat space to locate the 
black hole on the polar (positive z-) axis. The unit vector N normal to ^+ has Cartesian components A^* = 
(sin 6* cos sin 6* sin 0, cos 6'). The displacement from the black hole location to a point on has components 

r* = sin 6 cos (j>, i?+ sin 6 sin cj), i?+ cos 9 — D) = rn^ , 

where n is a unit vector. With q = D / R+ and x = cos 6*, the coordinate distance from the black hole to the point on 
J+ is 



r = R+ ^/l + q^ - 2qx. 
The scalar and vector products of the two unit vectors are 

1 — qx q sin 9 



N-n: 



N X n : 



(— sint/), coscj), 0) 



v/1 + g2 - 2qx g2 _ 2qx 

The direct contribution to the Bondi-Sachs mass aspect Ma from the conformal extrinsic curvature is 

3 T K 



(A/a)k = -j^Arr 



Substituting the expression in Eq. p2)) gives 



R 



(Ma)k = - 3 (N . n)^ - 1 -K[-] S • (N x n) (N • n) 



KC 
~6 

KR_ 
4~ 

K 

ml 



-—R^N'^AmN' 



R 



(Al) 



R4 



2 (N ■ P) (N • n) + (n • P) ( (N • n)^ - 1 



[2 (N • Q) (N • n) + (n • Q) (1 - 5 (N • n)') 
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The contribution of this to the Bondi-Sachs energy is the average of Eq. (|A1I) over sohd angle. Since the average of 
any term with an odd power of cos0 or sin0 vanishes, the spin term and the components of P and Q perpendicular 
to D do not contribute. We have 



(F ) 



KR^ 



— I dx- 



{l-qxY 



1 



{1 + q'^ - 2qxf'^ 2 7 {I + q2 - 2qx) 



dx 



\3/2 
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dx 



dx- 



X {1 — qx) 
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xil — qx) 



(1 + ^2 - 2qx) 



5/2 



dx 



— i dx- 



{x - g) (1 - qx) _ 
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dx- 



{x ~ q) 



2qx) 



{x - q) 



(1 + ^2 - 2qx) 



5/2 



— dx 



(x - g) (1 - qx) 
(1 + ^2 _ 2qx)^/2 



Evaluating the integrals with Mathematica, 



(^Bs)i 
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3-2g2 1 

4 2g-93 



4 



-q + O-O 
2 5q- 2q3 

3(1-92)2 ^3(1_^2)2 "3(1_^2)2 



which simplifies to 



(i?Bs)K = -^ + f (D-P) = |[-C + 2D.P]. 



(A2) 



The Bondi-Sachs momentum Pbs is the average over solid angle of Ma N . The portion of this directly due to the 
extrinsic curvature is best calculated separately for the components parallel and perpendicular to the displacement 
D of the black hole. The component parallel to D is obtained simply by inserting a factor of a; = cos 61 into each of 
the integrals in the above expression for {Ebs)k ■ Evaluating the integrals gives 



(^bs)k 



KC 
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KR, 
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15 15 
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which simplifies to 
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(A3) 



The perpendicular component has a contribution from the spin term in the extrinsic curvature. For instance, the 
x-component is 



(^'bs)k 



dx 



(1 - a;2) (1 ~qx) ^ I 
il + q^-2qx) 



3/2 



2 ''^^ 



(1 + g2 _ 2(;cc)^/' 2 



l-a;2 



K 

SR- 



KD^sy 

AR+ 



dx 



dx- 



{l - x"^) {I - qx) 1 
(\ + q^ -2qx)^'^ 2 

l- x^) [l-qx) 
{1 + q^ -2qx)'''^' 



dx 



'l-x' 



(l + g2_2gxf/2 



dx 



{l + q^ -2qx)^''^ 

[1 - x'') {1 - qxf 
(1 + g2 _ 2qxf''^ 



21 



Evaluating the integrals gives 



(-Pbs)k 
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px 
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Combine this result with Eq. (|A3[) to get in vector form 



(Pbs)k 



3 
K 



CD + P (1 - iDf ) +2(D ■ P)D + Q + CD x S 



(A4) 

The expressions in Eqs. (|A2[) and (|A4[) are only portions of the total Bondi-Sachs energy and momentum, respec- 
tively, so nothing can be concluded from these expressions about how the total Bondi-Sachs energy and momentum 
depend on the Bowen-York parameters. In particular, the total Bondi-Sachs energy is always positive even when the 
typically dominant first term in Eq. (|A2p is negative. Also, note that a displaced "Schwarzschild" black hole (with no 
boosts or spin) has, according to our numerical results, a net linear momentum in the opposite direction from that 
implied by the first term of Eq. (jA4p . 



Appendix B: Analytic calculation of total angular momentum 

For conformally flat initial data on CMC hypersurfaces, the expression in Eq. (j49p for the angular momentum at 
^+ is just an angular integral over the conformal extrinsic curvature A^'', independent of the conformal factor fl. 
The generalized Bowen-York expression for A^'^, as given in Eq. (|35|), is a simple sum over terms for each black hole, 
each of which contributes independently to the angular momentum. Therefore, we can consider each black hole and 
each term separately. For a given black hole, we can rotate the Cartesian coordinates in the conformal flat space to 
locate the black hole with some displacement D from the origin to put it on the positive z-axis. The unit vector N 
normal to this sphere at i? = the unit vector n directed away from the black hole, and the coordinate distance r 
from the black hole to a point on arc given in Appendix \X\ Also note Eq. (A3) for the scalar and vector products 
of N and n. Define q = D/R+ and x = cos6. 

The contribution from the first term in Eq. ([^^ is 



8^ 



i?4 



[3 (Nxn) (n-N) - (N X N)], 



which is identically zero since the integrand has only odd powers of sin (jj and cos < 
The "spin" term contribution to the total angular momentum is 
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(Bl) 



Consider separately the components of J along and perpendicular to the displacement of the black hole. The compo- 
nent parallel to D is 



dx 



(1 - qxY 



(1 + g2 _ 2qxf^^ 



2 3 - 2g2 2 1 



3 1 - 31 -q^ 



dx 
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(B2) 



These integrals (and those below) were done in Mathematica. The average over (f> kills the contributions of and 

sy to j§. 
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In the expression for the perpendicular component J^, there are non-trivial contributions from all three terms in 
the average over 0, with only terms proportional to sm'viving, 



J: 
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(B3) 



Similarly, = . The spin is unmodified by any factors depending on the ratio q = D/R^. 
Now consider first the normal boost terms in Eq. (|42p arising from the boost vector P, 
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sm 



dx 



(1 + 92 -2(7.t) 
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(B4) 



Again, the average eliminates all terms containing odd powers of sin (j) or cos (/). We see immediately that only 
contributes to Jf, and only contributes to and Jp. In particular. 
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(B5) 



and similarly for Jp. 

The other boost vector, Q, potentially could contribute to the angular momentum, because even though its contri- 
bution to the extrinsic curvature falls off more rapidly with coordinate distance from the black hole, is at a finite 
coordinate radius. The terms have a similar form to those in Jp, with 
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so in fact, the Q boost does not contribute. 

Our final result for the angular momentum of a generalized Bowen-York black hole at an arbitrary displacement D 
from the origin of a CMC hypersurface is 



J = S + D X P. 

Again, there are independent contributions from each black hole in the case of multiple black holes. 



(B7) 
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